MASS DETERMINATION WITH GRAVITATIONAL 

MICROLENSING 



Shude Mao^ and Bohdan Paczynski^ 
^ Max-Planck- Institute fiir Astrophysik, Karl-Schwarzschild-Strasse 1, 85740 Garching, 

Germany 

'\0 \ E-mail: smao@mpa-garcliing.mpg.de 

^Princeton University Observatory, Princeton, NJ 08544-1001, USA 
13-mail: bp@astro.princeton.edu 

Ph. 

. Subject headings: dark matter - gravitational lensing 

ABSTRACT 

> 
(N 

O . 

O . We present a simple toy model of the distribution of objects responsible for 
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gravitational microlensing. We use Monte Carlo simulations to demonstrate how 



\0 . difficult it is to determine the parameters of the lens mass distribution on the 

^ . basis of the observed distribution of event time scales. A robust determination 

P_|. requires ~ 100 events, or more, even if the geometry of lens distribution, and 



the lens kinematics are known. 



1. Introduction 

•i-H 

X. 

^ , One of the main objectives of the searches for microlensing events (Paczyhski 1996, 

and references therein) is to determine the mass function of the lensing objects. The 
determination of typical lens masses, or even their distribution function, were attempted by 
Alcock et al. (1993), Udalski et al. (1994), Zhao et al. (1995, 1996), Han & Gould (1996), 
and many others. 

We know that some lenses must be ordinary stars, but some may be brown dwarfs, 
and perhaps also planetary mass objects, and even more exotic things like black holes. 
It is recognized that even if all lensing objects were of the same mass M, they would 
give rise to a broad range of time scales to of the observed microlensing events. Although 
this is known to a relatively small number of people in the field this fact is not generally 
appreciated. Considering the broad interest in the nature of dark matter, and the relevance 
of microlensing searches to this problem, we feel it is justified to present a very simple 
model of the lens distribution and kinematics which allows some results to be obtained 
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analytically, and the rest can be calculated with one-dimensional numerical integrals. This 

model retains some of the most important generic characteristics: a very broad distribiition 
of to with power law tails towards very short and very long time scales for a delta function 
distribution of lens masses. Using Monte Carlo simulations we demonstrate how difficult 
it is to obtain accurate information about the mass function even within the framework 
of our simple model. In reality one is not sure what is the correct space distribution and 
kinematics of lensing masses, making the interpretation even more difficult. 

The model is described in the next section. The results of the Monte Carlo simulations 
arc presented in section 3, and the discussion of the results is given in the last section. The 
technical details of the model can be found in the appendix. 



2. A Model 

We consider a simple case that may approximate a real astrophysical system: the 
lenses, each of the same mass M, have their distances Dd distributed uniformly and 
randomly between the observer and a stationary source, which is at a distance Dg. The 
time scale of a microlensing event is given by 

to = ^, Vt = Vsmi, (la) 

where V is the relative spatial velocity between the lens and the observer, i is the angle 
between the velocity vector and the line of sight, and Re is the Einstein ring radius: 

C Us 

More technical details about our model are given in the appendix. In this section we present 
only the most essential information. 

The lenses in our model have a three-dimensional Gaussian velocity distribution with 
a one-dimensional dispersion of uy. A characteristic time scale for a microlensing event is 
defined as 

The probability distribution for event time scales has power law tails: 
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where 

i = ioAa, (3c) 

and the subscript S indicates that all the lenses have the same, i.e., we have a 5 mass 
function. The full probability distribution of event time scales can be calculated numerically 
and it can be approximated within 6% with a formula 

Psilogt) dlogt ^ X 10-l-44{logt+0.293)exp[-((logt+0.28)/0.304)2] ^^^g^^ 

where 

25.5 

c= — In 10 18.71, 

The first two moments and the standard deviation of this distribution are 

roo 

< logt >5= / psit) logt dt = -0.242, (5a) 
Jo 

roo 

< {\ogtf >s= / ps{t) {logtf dt = 0.116. (56) 
Jo 

as = [< (log tf>s-{<logt >sf] = 0.240. (5c) 

Naturally, formula (4) applies only to our simple model of the lens distribution and 
kinematics. However, the presence of power law tails is common to almost all distributions 
ever proposed (e.g., those in Zhao et al. 1996) and can be understood intuitively. The 
very long events are produced by the lenses which move almost along the line of sight. For 
these events the time scale t is inversely proportional to sini, i.e., t oc (cf. eq. la). 
The number of these events, A^, is proportional to the solid angle (1 — cosi) divided by t, 
i.e., oc (1 — cosi)/t oc oc t~^. On the other hand, the short events are produced by 
the lenses which are very close to either the source or the observer. For the latter case, 
t oc Re oc I^y^, and the number of these events is oc R^D^/t oc t^. The same argument 
applies when the lenses are close to the source. 

It is not likely that all the lenses have identical masses, and it is convenient to consider 
a power law mass distribution 

piM)di — )=A{ — ) di — ), for M^i„<M<M^ax, (6a) 



MoJ \MqJ \Mq 
where 

Mo = {M^r.M^^fl\ (66) 
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and we denote the logarithmic width of the mass function by 

(3 = log {M^^/M^^). (6c) 

The exponent a = —1.5 corresponds to an equal rate of microlensing events per decade 
of lens masses, while the case a = —2 corresponds to an equal contribution to the optical 
depth per decade of lens masses. For —2 < o; < —1.5, the optical depth (total mass) is 
dominated by the massive objects, while the event rate is dominated by the low mass 
objects. For a given mass range the standard deviation of logt is the largest for a — —1.5. 
For a <^ —1.5 the mass spectrum is dominated by low mass objects, and the event time 
scales have a small effective range, with the majority of events caused by lenses with 
M Mmin- In the opposite case, with a ^ —1.5, the microlensing events are dominated by 
lenses with M ^ Mmax- Therefore for a > —1.5, the standard deviation in logt approaches 
some asymptotic value when Mmin — > 0, and similarly for a < —1.5 when Afmax cxo. After 
some algebra one may show that the variance in the observed distribution of event time 
scales can be expressed as 

<t = ^l + ^, for a = -1.5, (7) 
and the asymptotic value is 

<t,oo = ^l + [2 (« + l-5) lnlO]-% for a ^ -1.5, /3 ^ oo, (8) 
where as = 0.240, as given by eq. (5c). 

The probability distribution functions for the logarithm of event time scale for any 
lens mass function can be calculated by convolving the mass function with the probability 
distribution for a 6 mass function as given in eq. (4). The probability distributions for 
(a, 13) = (-1.5,1), (-1.5,2), (-2.5,1), (-2.5,2), are shown in Fig. 1. AU time scales 
are in units of tcr corresponding to M = Mq = (MminMmax)^^^ (cf. eqs. 2 and 6). Also 
shown is the case for P — 0, i.e., with all the lenses having identical masses (cf. eq. 4). 
The shapes for mass functions with a narrow width (/3 ^ 1) are quite similar. This can be 
easily understood because the time scale of microlensing to oc M^/^, therefore one decade 
of mass range translates to a factor of three in the time scale, hence the shapes of these 
distributions still resemble that for a S mass function. 

3. Monte Carlo Simulations 

We present two types of Monte Carlo simulations. In both cases we adopt the model 
of lens distribution and kinematics as developed in the previous section. First, we consider 
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a perfect detection system, which can detect every event of any time scale, and a lens 
population with a moderately broad power law mass function. Next wc consider very broad 
lens mass functions and a detection system which is like that of the MACHO and OGLE 
collaborations, i.e., it has a significant detection efficiency over a relatively narrow range of 
event time scales. In both cases we make a large number of the Monte Carlo simulations 
to determine the accuracy with which the lens mass function parameters can be estimated 
with a moderately large number of detected events: n = 10, n = 100, and n = 1, 000. 

Wc tried both the moment method and the maximum likelihood method to estimate 
the parameters a and /5. Although the moment analysis is faster and in many cases its 
results agree with that from the maximum likelihood, it is less accurate and in some cases 
does not work at all. One of the reasons is very easy to understand. Imagine we have a 
random sample of microlensing events with variance smaller than — 0.240, as given by 
eq. (5c). Such cases are unavoidable in a finite sample drawn from a distribution with 
a > a^. The moment analysis fails for these cases. This reminds us of a similar old problem: 
what to do with negative trigonometric parallaxes of distant stars, the negative numbers 
being due to the unavoidable measurement errors? For this reason, we only present the 
results from the maximum likelihood method with a prior of equal probability in a and 
However, general discussions are given using the more intuitive moment analysis. 

3.1. Finite Mass Range, Perfect Detection System 

We adopt a power law lens mass function described with two dimensionless parameters: 
the power law exponent a, and the range of masses /5 = log (Mmax/^min) (cf- eq. 6). In 
addition, we have two scaling parameters: A and Mq = (MminMmax)^^^ (cf. eq. 6). Other 
things being equal the optical depth and the rate of microlensing events are proportional 
to A, while the average time scale of the events is proportional to . Out of the four 
parameters, the scaling parameter A can be fixed by the observed event rate, after which 
we are still left with three parameters. We can simplify the problem further by identifying 
a mass scale in the problem. In most cases, we assume such a scale is provided by a fixed 
upper mass limit, say Mmax = l-^o, but we also explore another possibility (see below). 

We assume that the lens detection system is perfect, i.e., it is 100% efficient for all 
event time scales. A large number of Monte Carlo samples are generated with the average 
number of events n = 100 and n — 1, 000. The actual number of events was chosen from a 
Poisson distribution, followed by a Monte Carlo selection of event time scales according to 
the model probability distributions, as shown in Fig. 1. 
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We selected two models with {a,P) — (—1.5, 1), (—1.5,2), each for two sample sizes: 
n — 100 and n = 1000. The la and 2a contours are shown in Fig. 2. It is clear that 
for the case of {a, (3) = (—1.5, 1), the power law index a cannot be determined very well 
with 100 events, the case for {a, (3) = (—1.5,2), however, is improved. The large difference 
between [3=1 and /5 = 2 is because for (3=1, the shapes, and therefore the variances, of 
the duration distributions for different power law indices are similar (cf. Fig. 1). The main 
difference between the models is the average duration, i.e., the first moment of the duration 
distribution. However, such difference can be accommodated by slightly changing the lower 
mass limit. The power law index a is then mostly constrained by the more uncertain higher 
moments, and hence the weak limit. 

The case with a broader mass function, = 2, is better constrained than that for /3 = 1 
with a sample oi n = 100 microlensing events, but even this case is very poorly constrained 
when n = 10, as shown with the dotted contours in Fig. 2. Note that the total number 
of microlensing events detected so far towards the LMC is smaller than 10 (Bennett et al. 
1996) and the spatial distribution and kinematics of the lenses are not known, making the 
case even more ill-constrained. 

The contours shown in Fig. 2 are based the assumption that the upper mass limit is 
known. Such a fixed upper limit sets the mass scale in the maximum likelihood method and 
translates the observed average duration (i.e., the first moment of the duration distribution) 
into a constraint on the lower mass limit. In reality, we do not know the upper mass 
limit well. In fact, this upper limit should also be found from the maximum likelihood. 
However, it is too time consuming to conduct such a three dimensional study. Instead we 
adopt a different approach: for a given set of {a, /?), we fix the mass scale such that it gives 
identical first moment in the duration distribution as a generated Monte Carlo sample. 
This reduces the three-dimensional problem into a two dimensional one, and should be a 
reasonable approximation. As the first moment has been utilized to set the mass scale, the 
determination of a and /3 is then effectively based on the second and higher moments, and 
these are subject to larger random uncertainties. The likelihood contours shown in Fig. 2 
confirms this for a = —1.5,j3 = 2 with n = 100 (dashed lines). Just as expected, the limit is 
relaxed - the allowed parameter space is increased by a factor of ^ 2. This is a consequence 
of the determination made with the maximum likelihood (or any other) method effectively 
using the higher moments of the event distribution. 

With the Monte Carlo samples, we can also estimate the optical depth of microlensing, 
r, which is proportional to the number of events observed per unit time per unit 
number of stars. The error of the optical depth estimate is inversely proportional to 
the square root of the number of events, with a coefficient somewhat larger 
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than unity, as the relatively few longest events are most important. This coefficient 
depends on the values of a and /9. For example, the coefficients are 1.3, 1.5, 1.3, 1.6 for 
(a,/?) = (—1.5, 1), (—1.5,2), (—2, 1), (—2,2), respectively. This is in good agreement with 
Han & Gould (1995). If we have 100 microlensing events, and if we can trust our model of 
the lens distribution and kinematics, then we can estimate the value of r with an accuracy 
of ^ 15%. 



Now we consider a very broad mass function with /? ^ 1, so the time scales of 
microlensing events cover the range A log to ~ 0.5/5. We adopt the efficiency of event 
detectability similar to those published by the OGLE (Udalski et al. 1994) and MACHO 
(Alcock et al. 1995a) collaborations. We approximate the detection efficiency with a 
formula: 



of detection sensitivity. The detectability window as described by eq. (10) has a width of 
about one decade at half maximum. Using eq. (2) we may relate this time scale to a mass 
scale 



Again we adopt a power law mass function as given by eq. (6). It is convenient to 
express the lower and upper mass limits in units of Mp with (5^-,^ = log (Mmin/Mp) and 
/^max = log (Mmax/^p)- The probability distributions of event time scales for very broad 
mass functions, with P — 10 and slopes a — —2, — 1.5, — 1, are shown in Fig. 3. Also 
shown is the efficiency of event detection as described with eq. (10). 

If the probability function of event time scales is much broader than the width of the 
detectability window as given by eq. (10) then most events cannot be detected. We shall 
consider two cases. First, a mass function extending very far towards low masses, i.e., with 
/3min <^ —1, and the upper mass limit somewhat larger than Mp. Second, a mass function 
extending very far towards large masses, i.e., with /3max 1, and the lower mass limit 
somewhat smaller than Mp. We shall investigate the ability of the detection system to 
measure /?max in the first case, and /Smin in the second case. 

In all the simulations we adopt the average number of detected events to be either 
n — 100 or n = 1, 000, the actual number being drawn from a Poisson distribution. 



3.2. Very Broad Mciss Range, Imperfect Detection System 





(11) 
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Assuming that /^min <S — 1 we can search for the remaining two dimensionless parameters: a 
and /3max using the maximum hkehhood method, with A being fixed by the event rate. The 
results of model calculations are shown in the right panel of Fig. 4 for (ct, /?max) = (—1-5, 1), 
and (—1.5,2). With 100 events, a is well constrained in both cases. The parameter /Smax, 
however, becomes not so well-determined for /3max = 2. 

Assuming that /9max S> 1 we can search for the remaining two dimensionless parameters: 
a and /5min using the maximum likelihood method. The results of model calculations are 
shown in the left panel of Fig. 4 for (a,/3min) = (—1-5, —3), (—1.5, —2), and (—1.5, —1). 
The parameter a is again well constrained for n = 100. In comparison, the parameter firava. 
becomes ill-determined when it approaches —3. 

Note that the models with /3min = —3 and /9max = 2 have contours slightly shifted with 
respect to the true model values, as marked with crosses. If the distribution of event time 
scales is approaching the effective edge of the detectability window then some finite samples 
of events generate second (and other) moments which are incompatible with the model. 
For such cases a method using the moments simply fails. The maximum likelihood method 
can better handle such cases, but has some trouble as well for n — 100, as indicated by the 
small shifts of contours in Fig. 4. This can be remedied with a larger sample size such as 
n — 1000 (thick dashed lines) in both cases. 

Given the determination of a and either /Smin or /^max, as well as the 'observed' number 
of events 'n' we may estimate the parameter A in eq. (6). The error of the estimate can be 
written as 

(7A ~ a (12) 

where a a is the standard deviation in the ^sample /^modei ratio. The coefficient a depends 
somewhat on a and /3, but is generally in the range between 1 and 1.6. 



4. Discussion 

The Monte Carlo simulations presented in this paper demonstrate that even if we 
have a full knowledge of the space distribution and kinematics of the lensing objects, or 
equivalently, if we know the relation between the lens mass and the distribution of event 
time scales (see eq. 4 for an example), it may still be difficult to obtain the parameters of 
the lens mass function. The problem is particularly serious when the width of the mass 
function is narrow (/9 ^ 1). In addition, if the mass function power law is quite different 
from —1.5, then it will be difficult to probe the high mass end when a <C —1.5, and the low 
mass end when a 3> —1.5 (cf. eq. 8). 
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Recently, Han & Gould (1996) obtained a rather tight hmit on the mass function of 
the Galactic disk using about 50 microlensing events toward the Galactic bulge. In their 
analysis, the upper mass limit is fixed at IOMq in a rather ad hoc manner. As we have 
shown in Fig. 2, such a fixed upper mass limit will make the determination of the other 
parameters appear more accurate. A more realistic approach is to estimate the upper mass 
limit using the maximum likelihood in order to avoid carrying our a priori assumption into a 
limit on the physical parameters. In general, if at least three separate parameters are to be 
estimated, for example Mq, a, and (3, then any method must effectively use the information 
contained in at least the first three moments. However, the higher the moment the more 
uncertain is its value as estimated from a relatively small sample of events. Therefore, the 
accuracy of the determination of any distribution parameter is lower in the 3-parameter 
determination than in a 2-parameter determination. Unfortunately, the smaller errors in a 
2-parametcr determination could be misleading, unless we have an independent and reliable 
information about the value of the third parameter. 

In our second example, very broad mass functions are sampled through a realistic 
dctcctability window (cf. Fig. 4). The current published sensitivity windows allow one to 
probe about three to four decades of mass range and we find that the "amplitude" A and 
the local slope of the lens mass function a can be determined reasonably well with ~ 100 
events. However, the estimate of the range of masses within the detectability window, /^min 
or Pmax, may be very uncertain with as many as n = 100 events, as shown in Fig. 4. 

If the mass function is broader than the dctcctability window then the total event rate 
or the total optical depth (or total mass) cannot be measured as the majority of events 
are outside the detectability window (Fig. 3), having either too short or too long to. This 
points to the necessity of broadening the detectability window, a task easy to accomplish in 
the near future. Some first attempts have already been made by the EROS (Ansari et al. 
1995) and MACHO collaborations (Bennett et al. 1996). 

How can we be sure that we have reached the low and the high mass end of the lens 
distribution? In principle this is simple: we have to broaden the detectability window and 
we should detect so many events that the generic power law tails in the distribution (cf. eq. 
3) become apparent. This may call for well over 100 events. Such a high number of events 
is within reach for the Galactic Bulge, where the rate appears to be very high (Udalski et 
al. 1994, Alcock et al. 1995b, c). The determinations of any parameters of the distribution 
of lens masses based on as few as 10 events are subject to very large uncertainties, as shown 
with the dotted contours in Fig. 2. 

In reality we do not know what the geometry of the lens distribution is, and what is 
the lens kinematics. While looking towards the Galactic Bulge the majority of lenses may 
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be in the Bulge itself (Kiraga & Paczynski 1994, Udalski et al. 1994, Paczynski et al. 1994, 
Zhao et al. 1995, 1996), or they may be in the disk (Alcock et al. 1995b, c). While looking 
towards the LMC the majority of lenses may be in our galaxy (Alcock et al. 1995a) or in 
the LMC (Sahu 1994; Wu 1994). Depending on the location and kinematics the relation 
between the lens masses and the event time scale to may be very different, and this leads 
to an additional major uncertainty on top of purely statistical uncertainty discussed in 
this paper. It seems rather difficult to disentangle all parameters one needs to describe 
the distribution and kinematics of the lenses and their mass function on the basis of the 
observed distribution of event time scale. 

Fortunately, the future studies of the variation of the optical depth with the location 
in the sky will help to identify the lens location. If lensing of the Galactic Bulge stars is 
dominated by the Bulge lenses then the optical depth should vary rapidly with the galactic 
longitude (Kiraga & Paczynski 1994). If such a variation is not present this will indicate 
that the lenses are located mostly in the galactic disk. If lensing towards the LMC is 
dominated by the LMC lenses then the optical depth should increase strongly towards the 
LMC center (Sahu 1994). If the optical depth is uniform over the LMC then the lenses 
must be mostly in our galaxy. This is very simple in principle, but it will take hundreds of 
events to establish beyond reasonable doubt. Such a large number will be readily collected 
towards the Galactic Bulge, but it will take a very long time for the LMC as the rate in 
that direction is very low (Alcock et al. 1995a). Once the geometry of the lens distribution 
is established it will be possible to develop reliable models of their kinematics, and infer a 
reasonable statistical relation between the lens masses and the event time scales to- 

This project was supported by the NSF grant AST93-13620 and by 
the "Sonderforschungsbereich 375-95 fiir Astro- Teilchenphysik" der Deutschen 
Forschungsgemeinschaft. We are very grateful to Dr. Peter Schneider for a critical 
reading of the manuscript. 
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A. Appendix 

We begin with the simplest lens model: all lensing objects have the same mass M, 
the same three-dimensional velocity V, and their velocity vector directions have a random 
isotropic distribution. The source located at the distance Dg is stationary, and the number 
density of lensing objects is uniform between the observer and the source. 

If all the events had identical time scales, then the number of microlensing events 
expected in a time interval At would be given as 

2 At 

N — — UgT — , for to = const, (^ll) 

TT to 

where 2/7r is the ratio of Einstein ring diameter to its area, in dimensionless units, Ug is 
the number of sources monitored, and r is the optical depth of microlensing. In fact there 
is a broad distribution of event time scales as lenses have transverse velocities in the range 
of < Vt < y, and distances in the range < Dd < Dg. A straightforward but tedious 
algebra leads to the equation 

Stt At f°° 
N^ — UsT — ^N / po{to) dto, {A2a) 



16 tm Jo 

where 



V J Dd=0.5Ds \ J V 

is the time scale for a microlensing event due to a lens located half way between the 
source and the observer, and moving with the transverse velocity V, and the probability 
distribution of event time scales is given by 



Po{to) dto = PT{T)dT : 



where 



-6T{1 + T^) + (3 + 2T^ + 3T^) In 



1 + T 



1-T 



dT, {A3a) 



T=—, < < oo. {A3b) 

The probability distribution diverges logarithmically at to — tm, however, the integrated 
probability converges. Also note that the probability density distribution satisfies: 
Po{l/T) = T^pq{T). Therefore, the probability that the event time scale is longer than 
T(> 1) is equal to the probability that the time scale is shorter than T~^. 

The asymptotic behaviors of the integrated probability are given by 

i^o(< T) = for T « 1, {A4a) 
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The power law tails are generic to almost all lens distributions ever proposed. Also note 
that only the first two moments of the distribution are finite: (to) and (tg), while the third 
and higher moments diverge. Therefore, it is convenient to use a logarithmic probability 
distribution defined as 

Po(logto) dlogto = (In 10) to Po(^o) dlogto, {A5) 

as all moments of this distribution are finite. 

Now we shall convolve the results obtained above with a three-dimensional Gaussian 
distribution adopted for the lens velocities: 



p{V) dV^(-) exp ( — , {A6) 



where (Ty is the one-dimensional velocity dispersion, related to the three-dimensional 
velocity dispersion by ay = Vrmsj V^- The total number of events expected from the lenses 
with such a distribution of velocities is given by 

A^.-y n.r-, [Ala) 

where 

ay \ ) 

and the subscript 5 indicates that all the lenses have identical masses, i.e., the mass function 
is a 5 function. 

After some algebra, the probability distribution of event time scales can be written as 
Ps{t) = o / cxp(-— ) privt) dv, t^—, — , (AS) 

Z Jo Z Gy 

where pt was defined with eq. ( A3a) . The asymptotic behaviors are given by 

211/2 

^^(> ^) = ^ ^' ^ » 1- (^9^) 

Now let us further adopt a general power law distribution of the lensing masses: 

p(M) dM oc dM, for M^in < M < M„,a^. (AlO) 
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The normalized mass function can be rewritten in a logarithmic form: 

p(logm) dlogm = (^^] m"+^ dlogm, m = Mq = {M^inM^^y/'^ , (All) 

\C{a) J Mo 



where 

pmax f In(mmaxMmin), if « = -1 . . 

^ ^ I fmit?-m^+„")/(l + a), otherwise ^ ^ 



mm V max 



After convolving with the mass function, the total number of events is given by: 

^=(32) "'^i^^iW' 

where r is the total optical depth by all the lenses with different masses and velocities, and 
the duration is normalized to 

= t.{M = Mo) = ^ {GM^\ (^,4) 
The event duration distribution is given by 

1 /-mmax ^f, 

P(^)- ^. / m>5(m-i/2^) dm, t = (A15) 

The asymptotic behaviors are given by 

211/2 _ Y\ 

011/2 I 2^ 1 

It is convenient to use a logarithmic probability distribution. A few examples of such 
distribution functions are shown in Figs. 1 and 3. 

De Riijula et al (1991) proposed to use the moment analysis to deduce the distribution 
of lens masses. After some algebra the first two moments can be expressed as 

< logt >=< logi >5 +^01, {Alia) 

< {logtf >=< (logi)^ >s < logt >s +^(/)2, {Al7b) 
<t = + ]i^2 - <^s = 0.240, {Al7c) 
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where the quantities with a subscript S were defined with eq. (5), and 

, ^m^" (log ^max) " ^Tl^^f +" (log m^j^) 1 

(1.5 + «)lnl0' 



^m^"(l0g^max)' " ^^f +" (log m^jn)' 1 ^/ . . 1 o, ^ 

'2 = 1.5+a _ M .^ + n,llnin '^1' " ^ ^^^^^^ 

max 



<'a^°-<?n+" (1.5 + a) In 10 

and 



4>l = 0.5 log (mmin^^max), « = -1.5, (Al9a) 

02 = ^ [(logmmin)^ + (I0gminin)(l0gmmax) + (log,niax )^] , « = -1-5. (Al%) 
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Fig. 1. — The probability distribution functions for the logarithm of event time scale are 
shown for four power law mass functions: (a,/?) = (—1.5, 1), (—1.5,2), (—2.5, 1), (—2.5,2) 
where a is the slope of the mass function and /3 = log (Mmax/Mnin). All the time scales are 
in units of corresponding to M = Mq = (Mmm-Z^max)^^^ (cf. eqs. 2 and 6). Also shown is 
the case for /3 = 0, which corresponds to all lenses having the same mass (cf. eq. 4). 
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Fig. 2. — The contours corresponding to one and two standard deviations are shown for 
different models of the mass function in the a — f3 plane, where a is the slope of the mass 
function and (3 = log (Mmax/^min)- The models correspond to {a,j3) = (—1.5, 1), (—1.5,2) 
as indicated by the crosses. The contours are for an average number of microlensing events 
n = 10 (dotted hnes), n — 100 (thin sohd lines), and n — 1,000 (thick solid hnes). The 
dotted and solid lines are for models with the mass scale set by fixing the upper mass limit, 
i.e., declaring it to be known. The dashed contours are for (q;,/3) = (—1.5, 1) with n — 100 
and with the mass scale set by the first moment (average) of the event durations. In all 
these simulations we assume that the observers can detect every microlensing event with 
the impact parameter smaller than its Einstein ring radius, i.e., that the detection system is 
perfect. 
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Fig. 3. — The logarithmic probabihty distributions for event time scales are shown for 
three very broad mass functions with (3 = 10 (/^min = — 6, /Smax = 4), and the slopes 
a = —2, — 1.5, — 1, respectively. The time scales are normalized to the duration at 
the peak sensitivity (cf. eq. 10). Also shown is the approximate shape of the efficiency curve 
of the MACHO and OGLE searches. Notice, that in this case the majority of lensing events 
are not detected as they fall outside of the relatively narrow detectability window. 



-19- 



1.3 



III I I I I 



1.4 



1.5 



-1.6 - 




1.7 



1.8 



_ ^ min 

" |S ■ 

1^ min 

1^ min 

^ min 

I I I I I I I I I I I I I I I I I 



1, n=100'' 

2, n=100 

3, n=100 
3, n=1000 




1 



---^max=l. n=100 

^max=2. n=100 

(3^«v=2, n=1000 

' max ' 

J I I I \ I I I I I I I I l_ 

1 2 

^max 



Fig. 4. — The contours corresponding to one and two standard deviations and an average 
number of microlensing events n — 100 and n — 1000 are shown in the a — (5 plane, where 
a — —1.5 is the slope of the model lens mass function. For clarity, only two sets of contours for 
n = 1000 are shown. The models on the left correspond to /3max = log {M^gj^/Mp) S> 1, i.e., 
the longest events are beyond the detectability window as given with eq. (10) and shown in 
Fig. 3. The contours are for the mass functions with /^jniii = log (Mmin/Mp) = —3, —2, — 1. 
The models on the right correspond to /3mm ^ —1, i-e., the shortest events are beyond the 
detectability window. The contours are for the mass functions with /3max = 1, 2. 



